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Using a Coulomb gas approach, we compute the generating function of the covariances of power traces for 
one-cut /3-ensembles of random matrices in the limit of large matrix size. This formula depends only on the 
support of the spectral density, and is therefore universal for a large class of models. This allows us to derive 
a closed-form expression for the limiting covariances of an arbitrary one-cut /3-ensemble. As particular cases 
of the main result we consider the classical /3-Gaussian, /3-Wishart and /3-Jacobi ensembles, for which we 
derive previously available results as well as new ones within a unified simple framework. We also discuss the 
connections between the problem of trace fluctuations for the Gaussian Unitary Ensemble and the enumeration of 
planar maps. 


I. INTRODUCTION 


Many textbooks on random matrices [4, 5, 47] begin with the derivation of Wigner’s semicircle law for the Gaussian Unitary 
Ensemble (GUE), 1 which is arguably the best known result in Random Matrix Theory (RMT). One of the most common proofs 
consists in introducing the moments GV, K = JV -1 'IVS'y, where 9 n denotes an TV x N GUE matrix and n is an integer, and in 
applying Wick’s formula to evaluate the averages (Gn,k)- In order to compute these expectation values, certain diagrams are 
associated to various contributions of Wick’s expansion; then, in the large N limit, only certain diagrams (the planar ones) survive, 
which can be counted explicitly. Eventually, one discovers that the limits limAr-^ ( Gn,k ) are non-zero only for k even and given 
by the Catalan numbers. This computation implies that the empirical density Pn(x) = A -1 JV <5( x — A;) of the eigenvalues {A,} 
converges in expectation to Wigner’s law p(x) = (27r) _1 x/4 — a; 2 . The simplicity of this result might lead to believe that the very 
same approach (known as the method of moments) should be easily extended to higher cumulants of Gn,k- Flowever, it turns out 
that this problem is highly nontrivial, with very few results available in the literature [11,15, 16, 25,26, 39,40,43]. To date, almost 
60 years after Wigner’s original paper [56], a self-contained explicit formula for the higher cumulants C v (Gat jKi ,..., G n,k,^ ) for 
generic v and {ft;} is still lacking 2 . The situation is even worse for general invariant ensembles, where the method of moments is 
not applicable. 

In this work we consider random self-adjoint matrices X ; .y whose eigenvalues { A } <E R. have the canonical distribution of a 
2D Coulomb gas at inverse temperature /3 > 0, namely 

i N 

dP^ttAfc}) = -- ln T(A i} ) TT dA ._ (1) 


Flere the single particle potential V(x) is assumed to be such that in the thermodynamic (large N) limit the Coulomb gas reaches 
a stable equilibrium distribution supported on a single bounded interval [a, b] of the real line. This class of one-cut ensembles 3 
includes the classical Gaussian, Wishart and Jacobi ensembles of RMT. For /3 = 1, 2 and 4, (1) is the joint law of the eigenvalues 
of real symmetric, complex Hermitian or quaternion self-dual invariant matrices, respectively. General ensembles parametrized 
by non-integer values of /3 > 0 can be realized from sparse real random matrices, as shown by Dumitriu and Edelman [24] and 
Killip and Nenciu [36]. Moreover, in view of certain applications to physics (see, e.g., [28]), it is worth studying the general joint 
probability density function (1) independently of a concrete matrix model realization. 

In this paper we consider the moments X^ K = 7V _1 TrX^ for matrices X^r with joint probability law of the eigenvalues given 
by (1), and we derive self-contained formulae for their covariances in the limit N —> oo. More precisely we study the quantities 


lim N 2 Cov(XN K ,XNe) = —a K e. 

iV—>• oo p 


( 2 ) 


The 1/^-dependence and the N 2 scaling of such covariances are customary in the fluctuations for 2D Coulomb gases [6]. The 
main result of the paper is a universal formula for the generating function of (2); it is universal in the sense that it depends only on 


1 For convenience and consistence with the formalism used for the problem of the enumeration of maps (see Appendix A), throughout this paper the GUE 
is defined by complex Hermitian N x N matrices Sjy whose diagonal and off-diagonal elements are normal random variables with the same variance, i.e. 
Sij = Xij , with Xij ~ Nc(0,1) for i < j , and xu ~ 3Mr(0, 1) for all i. 

2 C v (£, i,..., fu) denotes the cumulant of order v of the (not necessarily distinct) random variables £i,.. ., For instance C2(fi, £2) — (^1^2) — (£1) (£2)* 

3 The Cauchy-Stieltjes transform of the gas distribution has one connected cut on the real line of the complex plane, hence the name one-cut. 



2 


the support of the equilibrium density, but not on the potential V ( x ). This is a direct consequence of the macroscopic universality 
of the smoothed two-point function in RMT. Therefore, we can always write the covariances (2) in terms of a reference one-cut 
ensemble — the natural choice is the GUE. This allows us to prove a self-contained and manageable formula for the limiting 
covariances of an arbitrary one-cut ensemble (28). As particular examples, we derive known as well as new results for the 
Gaussian (Hermite), Wishart (Laguerre) and Jacobi /3-ensembles within a unified simple approach, unveiling the underlying 
structure. 

We will derive our results in the framework of the Coulomb gas analogy. The advantage is twofold. Firstly, the method is 
insensitive to the specific value of /3 > 0, so there is no need to produce separate proofs for each symmetry class. Secondly, 
it easily applies to matrix models whose underlying diagrammatic is cumbersome or unclear (such as the Jacobi ensemble). 
Moreover, the method is a good candidate for generalizations to higher order cumulants. 

Various methods have been employed to compute cumulants of moments of /3-ensembles. All of them are combinatorial in 
nature. There are several works on the Gaussian ensembles scattered throughout the literature of quantum field theory in the 
limit of large internal symmetry group [8, 12, 33, 34, 42]. Some of these results reproduce earlier findings by Tutte [51, 52] 
for planar enumeration of maps. Johansson [35] studied the global fluctuations of Hermitian random matrices and was the first 
one to prove a central limit theorem for /3-ensemble with polynomial potentials. More recently, in a series of paper devoted to 
the second-order freeness [16, 39, 40, 43], the relation between moments of some random matrix ensembles and non-crossing 
partitions has been employed to compute covariances to leading order in N on complex and real Gaussian and Wishart matrices. 
Dumitriu and Edelman [25] and Dumitriu and Paquette [26] studied the fluctuations around the equilibrium density in Gaussian, 
Wishart and Jacobi ensembles at generic Dyson index /3. Their approach was based on the realization of these ensembles in 
terms of sparse matrices. The computation of cumulants of moments for some /3-ensembles is possible in principle by solving a 
recurrence relation associated to the hierarchy of loop equations [10, 11, 15]. 

The plan of this paper is as follows. We start with some preliminaries and a summary of the main ideas of the Coulomb gas 
approach (Section II). In Section III we present our main result. In Section IV we analyse the classical /3-ensembles: the Gaussian, 
Laguerre and Jacobi /3-ensembles. Finally, in Section V we discuss some open problems. The paper is complemented by two 
appendices. In Appendix A we review the connections between the problem of moments fluctuations, enumeration of planar maps 
and non-crossing pair partitions. Appendix B contains the derivation of a technical result. 


II. 2D COULOMB GAS AND CUMULANTS OF LINEAR STATISTICS 

Here we briefly summarize (without derivations) some aspects of the 2D Coulomb gas approach, which are more directly 
relevant to our work. The starting point of the method is the analogy with a gas of Coulomb charges pointed out by Wigner and 
Dyson [27, 56]. The joint law Pjv,/ 3 ({A/c}) in (1) is the Gibbs-Boltzmann measure of a 2D Coulomb gas confined in a interval 
Aciat inverse temperature 8 > 0 (the Dyson index) with a single-particle potential V (A) bounded from below and finite for 
As A (we adopt the usual physical convention that probabilities vanish in regions of infinite potential). For simplicity we assume 
V ( x ) convex with superlogarithmic growth at infinity. It is known [35] that under these assumptions on V(x) the Coulomb gas 
is supported in the thermodynamic limit on a bounded interval [a, b\. These assumptions on V ( x ) are purely technical, and in 
fact the correct condition should be the stability of the one-cut configuration under small analytic perturbations of the external 
potential (this condition goes under the name of ‘off-criticality regime’). 

In the thermodynamic (large N) limit it is convenient to describe the positions (Ai,..., Xjy) of the gas particles (the eigenvalues 
of the /3-ensemble) in terms of the empirical density of the Coulomb gas, 

1 N 

Pn{x ) = — ^2 5{x - A i) . (3) 

2=1 

The partition function of the gas may be written in terms of the empirical density in the following way. First, notice that for large 
N the energy of the Coulomb gas admits an integral representation 

-^ 1 n|A i -A,|+iV^U(A i ) = 

iyij i 

in terms of the mean-field energy functional £ defined as 

= da(x)da(y) In \x — y\ 

^ JJ x^y 

Using this energy functional we have 

r r N 

Zn,0 = / V\pn] / , 

^ J 2—1 


N 2 Z{ Pn ]+o(N 2 ) 


da(x) V ( x ) . 


(4) 


(5) 


( 6 ) 
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where the functional integration f D[p,y] is restricted to the probability measures compatible with the Coulomb gas configuration 
{ Ai}. A spectral linear statistics (or linear statistics for short) is defined by 

N 

A n = N~ l Y, /(Ai) = 

i =1 


da; p N (x)f(x), 


(7) 


where f(x) is a suitably chosen test function. Sometimes (7) is referred to as sum function on the particles positions {A;}. The 
average of Ayr is given by 



dx (pn(x)) f(x). 


( 8 ) 


In a similar way one computes the variance of A y, 

N- 1 Eh /t(A<) and = TV” 1 Eh hi^i) as 


or in general the covariance between two linear statistics 


A 


(i) __ 
N ~ 


Cov(A^ } , = J dxdyCov(p N (x)p N (y))f 1 (x)f 2 (y) ■ 


(9) 


The same reasoning can be extended to higher order cumulants. In general, the mixed cumulant C„ ( A $, ,..., ) of v 

linear statistics (not necessarily distinct) can be recovered by integrating A$ = /V -1 Eh against the cumulant of the 

gas density at v points C v (pn{x\), pn{x 2 ), • • •, Pn(x v )), also known as v-point connected correlation function: 


C (d 1 ) d 2 ) 


.4?) 



C v {pn{x\),Pn{x2), ■ ■ -,Pn(x v )). 


( 10 ) 


For v = 1 and v = 2 this corresponds to (8) and (9), respectively. In the asymptotic regime as N —> 00 , the connected correlation 
functions of the gas are generated by the free energy — log Zn,p by taking a suitable number of functional derivatives with respect 
to the external potential. More precisely, for every positive integer v: 


C v (p N (xi),p N (x 2 ), • • •, pn{x v )) 



n 


6 

8V(xi) 


In Zn,/3 [1 + o(l)]. 


( 11 ) 


For instance, the 1-point and 2-point connected correlation functions in the limit as N —t 00 are given by 




/31V 2 5V(x) 


Co v(p N (x)p N (y)) = 


1 


S 1 In 2 >N,p 


[l + o(l)] 


JN 2 J SV{x)SV(y) 

Under the assumption that the limit limy^oo (py(x)) = p(x) exists, from (11) we see that 


( 12 ) 

(13) 


lim 

N—fOO 


1V 2 + 1 ' , C v {p N {x 1 ), pn{x 2 ),..., pn(x v )) 


v—1 


"v—1 

n 


6 

8V(xi ) 


lim (p N {x v )) 

N—foo 


(14) 


exists too 4 . Hence, in principle a complete answer to the general problem of computing spectral linear statistics in the large N 
limit could be achieved if all functional derivatives with respect to the external confining potential were known. At present the 
complete solution to this problem is out of reach. Fortunately, in some circumstances we have a partial but explicit answer. In the 
following we focus to first and second order effects: 


P{x) = lim (pn(x)} 

N—foo 

1 So(x) 

X(x,y) = - A linr o iV 2 Cov(p Ar (x)p A r(2/)) = psY^j)' 


(15) 

(16) 


4 The last equality for v = 2 was first recognized by Beenakker [6] in the random matrix approach to quantum transport in chaotic cavities. 
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The quantity p(x ) is the density of states, or equilibrium density, while X(x, y) is known as smoothed two-point connected 
correlation kernel (see [7, 13] for details). Therefore, the limiting mean and covariances of smooth linear statistics (the lowest 
order cumulants) are 


lim 

N—fOO 


( a n) 


j d xp(x)f(x) 


(17) 


lim N 2 Cov(A^,A^) = — [ dxdyX(x,y)f 1 (x)f 2 (x). 

v —^OO J 


(18) 


From now on, we will restrict ourselves to one-cut ensembles, i.e. we will always assume that the potential V(x) is convex and 
with superlogarithmic growth at infinity. Under these hypotheses, the density of states p(x) (15) exists, is absolutely continuous 
and supported on a single bounded interval supp p = [a, 6], (a < b), and is the unique probability measure that minimizes the 
mean-field energy functional £ in (5). We notice that (3 does not appear in the minimization problem, and hence the density of 
states (and in particular its support) does not depend on (3. Remarkably, under the condition that the support is a single interval, 
p(x) and X(x, y) have the closed expressions 


p(x) 

X(x,y) 


*-a<x<b 


l 

n \J 0 - a)(b - x) 

1 1 a<x,y<b 


1 +~P 

7T J 


rb dy v '{y)V(y- (i )( b -y) 


y-x 


d sj(x — a)(b — x) 


/3tr 2 \/(y — a)(b — y) dx x-y 


(19) 

( 20 ) 


where P denotes Cauchy’s principal value. 

The first result (19) is usually referred to as Tricomi’s formula [49]. The beautiful identity (20) is a consequence of Tricomi’s 
formula and (16), and has been independently discovered by several authors [3, 7, 13, 14, 29]. It is worth emphasizing that 
while the one-point density (19) depends explicitly on the potential V(x), the smoothed two-point connected correlation kernel 
(20) enjoys a very surprising feature: it depends only on the edge points a, b of the spectral density. The dependence of X(x, y) 
on the details of the model only through the edges of the spectral density is often referred to as to long-range (or macroscopic) 
universality. Conversely, the microscopic universality emerges on a scale 0(l/iV) of the interior points of the limiting support of 
p(x) and is governed by the celebrated sine-kernel [13]. The universality of X(x, y) is a key observation that eventually leads to 
our main result, presented in next section. 


III. MAIN RESULT 


Let Xm, k — N 1 TrX'y denote the moments of a generic one-cut /3-ensemble of random matrices. For any fixed positive 
integers k and t, the limits 


exist and are given by 


a K = lim (X Nli ) 

N—foo 

\ot K y = lim N 2 Cav(X N K ,X N t ) 

p Af->oo 


a K ,t 


dx p{x)x K 

rb 


= —P 


dxdy, 


(x ~ a)(b - x) k « j 
(y-a)(b-y)y-x y 


( 21 ) 

( 22 ) 


(23) 

(24) 


respectively. The identities (23)-(24) are a specialization of (17)-(18)-(20). While the limiting average a K depends explicitly on 
the probability measure that defines the ensemble through V (x), and hence p(x) (given in (19)), the covariance (l//3)a K .« does 
not. Indeed, the following theorem is universal, as it depends only on the fact that the density of states is supported on a single 
bounded interval. 


Theorem. Let Xn belong to a one-cut (3-ensemble ((3 > 0) and let supp p = [a, b] denote the support of the density of states. 
Then, in the limit N —> oo the generating function of the covariances (22) is 


Fu 


1 


i < 


(,,0 = s ‘ = 


k ,£=0 


Hz- C) 5 


2 abzC^ — {cl -\- 6) (z H- £) H- 2 
2 V (1 - az)(l - bz)( 1 - aC)(l - bQ 


- 1 


(25) 
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This is the main result of this paper. 

Remark. It is worthwhile to emphasize the following properties of F[ a b ] (z, (). 

1. F[ a b ] ( z , () = F[ a b ] ((, z), hence a K ^ = a^ iK , which expresses the fact that the covariance is a symmetric functional. 

2. One can verify that liiru^o F[ a ,b] ( z > 0 = 0- Therefore o K j = 0 whenever k or i are zero. This is a consequence of the 
fact that 2tjv ; o = 1 is a constant. Moreover, F a j; (z. Q is continuous in (z, Q and analytic in each variable (keeping the 
other variable fixed) in a neighbourhood of the origin that does not intersect the cuts of the square roots in the denominator. 
If a function of complex variables is continuous and separately analytic in each variable, then it is also jointly analytic 
(Hartogs’ theorem [32]); therefore, Fj a h ] ( z , C) is analytic in a neighbourhood of (z, £) = (0, 0). 

3. The limiting covariances are proportionals to 1//3, which is true whenever V(x) is independent of 3 (see, e.g. [7, 28]). 

4. The generating function (25) satisfies the homogeneity identity (2, £) = Fi a u(tz,t£) for all real t. Therefore, 

under the map (a, b ) K»• (fa, tb) the limiting covariances transform as K > t K+e a Kt £. This implies that the correlation 

coefficients r K> ^ = ce K ,e/y/(XK.,K,ote,e are invariant in absolute value. More precisely, 

r K ,e (sign(f)) K+ V K) £ as (a,b) n> (ta,tb). (26) 


5. When the support of the equilibrium density is centred at the origin, i.e supp p = [—b, b] with b > 0, then the generating 
function simplifies considerably; 


F± b (z,Q 


1 zC 1 - b 2 z( 

P(z-C ) 2 ^/(l-b 2 z 2 )(l^¥C) 


(27) 


Furthermore, it satisfies F±b(z, Q = F±b(—z, —£); as a consequence a K} ( = 0 whenever k ^ i mod 2. 

6. The support [a, b] of the density of states can always be centred at the origin by a shift A i—>• A — m, where m = (a + b) /2. 
This corresponds to a translation in the matrix model >-> Xjv — rn/jv where 1^ is the identity matrix. It is known that 
the 2D Coulomb gas interaction is invariant under conformal transformations [17]. In other words, the covariance structure 
depends uniquely on the fact that the support is a single bounded interval. The translation Xjv >—>• X,v — ml m induces the 
change of variables x i—► x — m and y > y — m into the integral in the r.h.s. of (24), which leads to the relation 


= m 


k,-\-£ 


EE 

p—0 q —0 


PJ w 


\ P+Q 


— 
2m J 


orz 


(28) 


where to = (a + 6)/2, L = (b — a )/2 and ct| are the covariances of an equilibrium density supported on [—2, 2] (e.g. of 
the Gaussian ensemble). The are explicitly given in (42). Thus, (28) gives a tool to evaluate the limiting covariances 
for an arbitrary one-cut /?-ensemble. This is somehow analogous to the usual process of centring the variables in the theory 
of probability. 

7. If a K> e £ N for all k, £ then max(|a|, |6|) > 1. This is a necessary condition to have positive integers a K ^. 


In what follows we present the proof of (25). As in [19], the real integral (24) can be lifted to a double integral in the complex 
plane, 


a K( = lim lim — 

e->0+ e'->0+ 7 T z 




(z — a)(z — b) K 
(( ~ a)(( ~ b) ( - z 


(*z 


K—l 


(29) 


where T e , I V are two clockwise oriented contours enclosing the cut [a, b] at distance e and e' respectively (see Fig. 1). The first 
integration with respect to ( is performed with z fixed inside the cut. In this way the singularity (, = z becomes irrelevant. One 
verifies that the result of the first integration (a function of z) is analytic in a neighbourhood of the cut and hence integrable on the 
contour F f enclosing the cut. With this prescription in mind, for n.l > 0, we get using residues 


cuk ./ = k x Res 


( (z-a)(z~b) 1 t i 

\]J (C-a)(C~b)z-C 


;C 


= OO, 2 = 



(l — az)(l — bz) C 1 1 


(30) 


= k x Res 


(1 - <)(1 -bQz-C C^ +1 z K+1 


;C = 0,^ = 0 


(31) 
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lim 

e'— >-0+ 
z€.(a,b) 


ii) lim 

e->0+ 


IV 


rx 



f -*—• 


/ T 
e 

a z b 

(- 

Rc 


_1 

->—N 

v. 

e 

n 

a b 

(- 

? 

Re z 


FIG. 1: Contours of integration and scheme of the prescription used to evaluate (29). First the integral and the limit e —> 0 + in the ((-plane with 
2 inside the cut [a, b]. Then the integration in the 2 -plane. 


From (31) we have 




k x coefficient of z K f >: in the expansion of 


( I (1 - az)(l - bz) C \ 

\V _ “CX 1 


about (z, () 


= coefficient of z K ( £ in the expansion of 2 d z 


( / (l-a,)(l-^) C \ 

\\j (1 - 0(1 -b()z-(J 


about ( 2 , C) 


(0,0) 

(32) 

(0,0). 

(33) 


In the last line we have used the observation that if C(z) = c K 2 K , then the generating function of k c k is zd z C(z). Hence, the 
generating function F(z, <) = (1/X) J2 K ,e> 0 is given by 


F[a, 6] C) 


1 

X 


x analytic part at ( 2 , £) 


(0, 0) of the function 2 d z 


( / ( 1 -Q(i-fe) C \ 
\V X _ 0(1 - K) z- Cy 


(34) 


Carrying out the differentiation and isolating the analytic part, one gets the final expression (25). 

We also have an alternative derivation of (25), which relies on the chain of loop equations for the u-point resolvent 


W(z 1 ,.. .,z v ) = C v 


—Tr--- —Tr — 

N 2i-Xjv ’N z v 



(35) 


When the equilibrium density has square root singularities at the edges of the support (soft edges) these quantities admit a power 
expansion in 1 /N (we refer to [11, 15] for a more detailed exposition) and, under certain hypotheses, one can compute the ri-point 
resolvents recursively order by order by solving a hierarchy of loop equations (for an algorithmic approach see [10]). Within this 
formalism, the generating function (25) can be expressed as -Fj a .h]( 2 , C) = zC,w!f\\/z, 1/C) where W.^ is the leading order in 
1/iV of the 2-point resolvent. Using the notation in [10], we have 


F[a,b] (-^5 C) 

uj(x,y) 


2 


/?2C 


U) 


1 1 



(36) 

(37) 


where j a ^ is the Joukowski transformation defined as j a ,b(z) = + Xp { z + z ) (i n (36)-(37) j a l and j' a b denote the inverse 

map and the derivative of j a> b, respectively). 


IV. APPLICATIONS TO THE CLASSICAL /3-ENSEMBLES 
In this Section we focus on the classical /^-ensembles: the Gaussian, Wishart and Jacobi ^-ensembles. 

Fluctuations of traces in the Gaussian and Wishart /3-ensembles were considered by Dumitriu and Edelman [25] in the 
framework of the tridiagonal realization of /3-ensemble that they introduced [24] (for /3 = 1,2 the tridiagonalization of random 
matrices had already appeared earlier in [50] and [44]). Using a combinatorial technique simplified by the sparseness of their 
matrix models, they managed to obtain formulae for the averages and covariances of the moments (among other results). Their 
formulae are insensitive to the particular value of /3 > 0 (since /3 appears as just a parameter in the distribution of off-diagonal 
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terms). Later, Dumitriu and Paquette [26] investigated the very same statistics on the Jacobi ensemble. The matrix realization of 
this ensemble is not tridiagonal [36]; this turns out to be a serious obstruction to obtain explicit formulae. 

Here we show that our approach allows to treat these cases on the same footing: neither the specific matrix model nor the 
concrete matrix representation do matter. We will recover previous results, and produce new self-contained formulae for the 
limiting covariances of the Wishart and Jacobi /3-ensemble as particular cases of the general formula (25). 


A. Gaussian ensemble and symmetric intervals [a, 6] = [— L, L] 


The density of states of the /3-Gaussian ensemble is symmetric with respect to the origin. It is, therefore, a natural place to start. 
The eigenvalues of /3-Gaussian matrices Sw are distributed according to 


N 


dPiv.MUfc}) oc |A({Afc})I' 0 [] e-^^dAfc, 


(38) 


1 


where A({A fc }) = Ilj<fc (Afc — A j) is the Vandermonde determinant. Eq. (38) is easily recognized as the canonical measure of a 
2D Coulomb gas in a quadratic potential V(x) = x 2 /4. When /3 = 1,2,4 (38) is the joint probability density function of the 
eigenvalues of real symmetric (GOE), complex Hermitian (GUE) and quaternion self-dual (GSE) gaussian random matrices, 
respectively. 

For any /3 > 0 the average of the empirical density (3) of the eigenvalues {A;} converges in expectation to Wigner’s semicircle 
law. 


» w = 


|x|<2- 


Therefore, the averages of the moments a% = lim/v^oo ( Gn,k ) converge to the Catalan numbers 

for k = 0 even 
for k = 0 odd. 

The limiting density pg is supported on the interval [a, b] = [—2, 2]. Thus, formula (25) reads 

1-4 zC 


r 2 , 

( k \ 

k + 2 ' 

viJ 

lo 



F±2( z ,C) 


y/(l-4z*)(l-4?) 


- 1 


(39) 


(40) 


(41) 


The coefficients a/' , of this double series can be first extracted by taking fcth and (th order derivatives of (41) symbolically, 
and recognizing an easy pattern. This leads eventually to the well-supported ansatz 


a* e = lim pN 2 Cov{G N , K ,G N ,t) 

’ IV —Vrvi 





if k = £ mod 2 
if ac ^ £ mod 2. 


(42) 


where is the largest integer smaller than or equal to x. Then, it remains to prove that the guessed form (42) is indeed correct 
by resumming the series e ^Cf' - this is outlined in Appendix B. 

The first values of these covariances are given explicitly in Table I. We mention that (42) can be tracked back to a work by 
Tutte [52] in combinatorics (see Appendix A). 

More generally, using the homogeneity property of F[ 0 ^(z, C) ( see Remark III), for any one-cut /3 ensemble Xn whose 
limiting spectral density is supported on a single symmetric interval [— L, L] we have 


lim N 2 Cov(X N K ,X N e ) 

JV—foo 


1 

P 





For L = 1 this expression agrees with the formula derived by Dumitriu 
combinatorial proof. 5 


if [a,b] = [-L,L\. (43) 


if k = £ mod 2 
if mod 2 

and Edelman [25] with a much more complicated 


5 In [25] a different choice of the scaling of the eigenvalues leads to an equilibrium density supported on [—1,1], 



lim N 2 Cov (Gjv.k, Gn/) 

N —¥00 



- 2 

0 

6 

0 

20 

0 

70 

0 


0 

4 

0 

16 

0 

60 

0 

224 


6 

0 

24 

0 

90 

0 

336 

0 

1 

0 

16 

0 

72 

0 

288 

0 

1120 

p 

20 

0 

90 

0 

360 

0 

1400 

0 

0 

60 

0 

288 

0 

1200 

0 

4800 


70 

0 

336 

0 

1400 

0 

5600 

0 


0 

224 

0 

1120 

0 

4800 

0 

19600 


TABLE I: A block of limiting covariance structure (1 < k, l < 8) for the /3-Gaussian ensemble from (42). 


It is worth emphasizing that all the ensembles with a symmetric support of the equilibrium density share the same L- and 
/3-independent correlation matrix (see Remark III), 


r K ,e = 




2 -- if k = £ mod 2 n r t n 

K + e if [a,b\ = [-L,L\. 


V a ^ a W l 0 if k ^£ mod 2 


(44) 


B. Wishart ensemble and intervals [a, b] = [0, 2L] 

The same direct method that we used for the Gaussian /3-ensemble can be applied to the /3-Wishart (Laguerre) ensemble of 
positive semi-definite random matrices W/v > 0, characterized by a joint probability density of the eigenvalues 

N 

dP^({A fc })a |A({A fc })|^ JJ Af (c_1)JV e-^ JVAfc 0(A fc )dA fc , (c > 1). (45) 

k =1 


For /3 = 1,2,4, this ensemble can be realized as Wjv = ^9atxmS]vxM’ where S NxM denotes a Gaussian N x M matrix with 
i.i.d. real, complex or quaternion entries, respectively. The parameter c > 1 in (45) is the ratio M/N, which is assumed to remain 
constant as M, N — > oo. We notice that (45) is the canonical measure of a 2D Coulomb gas in the (non-polynomial) external 
potential V ( x ) = x/ 2 — (c — l)lna: for nonnegative x and infinite on the negative half-line. 

The average of the empirical spectral density of this ensemble converges to the Marcenko-Pastur law with parameter c > 1: 

Pw( x ) = ^\/(x-a)(b-x) 1 z 6 (o,f,), with [a, 6] = [(1 - \/c) 2 , (1 + v/c) 2 ] . (46) 

The limit on average a // = linrv^oo {Wn,k) of the moments I-'F/v.k = N~ 1 TvW] v exists and is described in terms of the 
Narayana numbers: 


= X Nar(«;,p), Nar(/c,p) = ^ ^ K _ ^ ^ ^ . 

Formula (28) combined with (42) gives 



o <q<e 
p—q mod 2 




(47) 


(48) 


Few values of a^ are reported in Table II. 

The /3-Wishart ensemble was also investigated by Dumitru and Edelman [25] within the framework of the tridiagonal realization 
of the ensemble. Their proof requires an ad hoc combinatorial analysis, which is different from the /3-Gaussian ensemble. Their 
formula (see [25], Claim 3.15.2) is considerably more involved than (48), and besides, we found that it contains few misprints. 
When the support is the interval [0, 2 L\, the generating function of the a K y simplifies drastically 


1 z( 1 -L(z + () 

y/(l-2Lz)(l-2LQ 


F{0,2L](z,0 


(49) 
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lim N 2 Cov(W n , k , W N , t ) 

N—too 


1 

0 


2 c 

4 (c + c 2 ) 

6 (c + 3c 2 + c 3 ) 


4 (c + c 2 ) 6 (c + 3c 2 + c 3 ) 

4 (2c + 5c 2 + 2c 3 ) 12 (c + 5c 2 + 5c 3 + c 4 ) 

12 (c + 5c 2 + 5c 3 + c 4 ) 6 (3c + 24c 2 + 46c 3 + 24c 4 + 3c 5 ) 


TABLE II: Block of the limiting covariance structure (1 < k,£ < 3) of the /3-Wishart ensemble from (48). 


Extracting the Taylor coefficients of this formula involves lengthy calculations. However, for c = 1 the moments of Wishart 
matrices are related in a simple way to the moments of Gaussian matrices and therefore one obtains a self-contained formula for 
the covariances, 

Jim ^Cov(X„,,*„,) = 1 (^)"'^( 2 “; ‘)( 2 V ')■ *M = [».“]■ W) 

Setting 2L = 4 we get the covariance structure of the symmetric (c = 1) Wishart ensemble. 

The correlation coefficients are independent of the length of the interval and are given 


rt 

r „ p = —= 2 - 

’ n + £ 


if [ a , b] = [ 0 , 2 L\. 


(51) 


In Fig. 2 we report Monte Carlo simulations of correlation coefficients for ensembles of type [0, 2 L\. From (50) we can derive the 
large k. £ behaviour of the covariances 


lim A+Cov (X N ' K ,X N j) 

N^oo 


(2L) K+e Vkl 
f3ir k + £ 


if [a, b\ = [ 0 , 2 L\, 


(52) 


in agreement with an earlier prediction of two of the authors of this paper [17]. 


C. Jacobi ensemble 

As last example we consider the /-Jacobi ensemble. These matrices satisfy the constraint 0 < 3n < / and their eigenvalues 
0 < A fe < 1 have joint probability density function 


N 


cffV l/ 3 ({A fc }) a |A({A fc })| /? H A^ 71 (1 - X k ^ N ^9(X k (l - X k ))d\ k , ( 7l , 2 > 0). 


k =1 


The density of states of this ensemble is 


ftlO) = 2to(1 2 -X) ~ a ^ b ~ ^ 


It is compactly supported on the interval [a( 71 , 72 ), &( 7 i, 72 )], whose edges ; 


1(71,72) = 


V 72 + 1 - a /(71 + l)(7i + 72 +1) 

71+72 + 2 


H 71 , 72 ) = 


V 72 +1 + Vhi + 1 ) (71 + 72 +T( 

71 + 72 + 2 


(53) 


(54) 


(55) 


By inserting (55) into (25) one obtains the generating function of the limiting covariances of Jn,k = TrJ'/ for any value of 
71,2 > 0. The a 3 f’s can then be evaluated explicitly via (28) as 


r 7i + 7 i72 + 2(71 +72 + I) 
(71 + 72 + 2 ) 2 


a K ,t = 


-1 k,-\-£ 


E 

0<p<K, 

0 <q<£ 


\/ (7i + 1)(72 + l)(7i +72 + 1) 
7i + 7 i 72 + 2(71 +72 + I) 


p+q 


1 c 

pj W 


QZ. 


(56) 


where ajj g is given in (42). Furthermore, in the particular case 71 = 72 = 0, we have [a, 6 ] = [0,1]; hence, using formula (50) 
we arrive at 


a 3 t = lim (3N 2 Cov {Jn, k , Jn,i) = 4 1 

’ N—>oo 


k£ /2k - 1\ / 2t - 1 


K + £ V K 


+ 71,2 = 0 ([a, b] = [0,1]). (57) 
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FIG. 2: Limiting correlations r K ,e = a K ,r/ yJa. KtK .on t i of the moments for the Wishart ensemble W n with c = 1 and the Jacobi ensemble 3n 
with 7 ! = 72 = 0. The densities of states pw and pg are supported on [a, b] = [0, 4] and [a, 6] = [0,1], respectively. According to Remark III, 
the correlation coefficients r K 7 are the same for these ensembles. Here we show a comparison of the prediction (51) with numerical simulations. 
The size of the matrices Wjv and 3n is N = 50, the sample size is n = 10 4 , k = 1,..., 8 and t = 1,..., 5. Dotted lines connecting the 
theoretical values are depicted to guide the eye. 


Few values are reported in Table III. 

Dumitriu and Paquette (see [26], Eq. (23)) derived an integral formula for a 3 K that can only be evaluated numerically; one can 
then check that the numerical integration gives the same values as formula (56). Their proof, however, does not extend to the 
limiting values 712 = 0 , which are covered by ours. 


9 1 

lim N"Cav (Jn,k, Jn,i) = -7 

JV->-oo ’ (3 


1 

1 

15 

7 

105 

8 

8 

128 

64 

1024 

1 

9 

9 

35 

135 

8 

64 

64 

256 

1024 

15 

9 

75 

75 

4725 

128 

64 

512 

512 

32768 

7 

35 

75 

1225 

1225 

64 

256 

512 

8192 

8192 

105 

135 

4725 

1225 

19845 

1024 

1024 

32768 

8192 

131072 


TABLE III: A block of the limiting covariance structure (1 < k, l < 5) of the Jacobi ensemble for 71,2 = 0 from (57). 
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Y. CONCLUSIONS 

Using a Coulomb gas approach, we have derived a formula (25) for the generating function of the covariances of the moments 
in the thermodynamic limit N —> oo for one-cut /3-ensembles of random matrices. This result is universal, as it depends only on 
the support of the limiting spectral density, but not on the potential V ( x ) that defines the ensemble, nor on the inverse temperature 
/3; this is a direct consequence of the universality of the smoothed two-point spectral correlation function. Given a /3-ensemble 
whose equilibrium density has support [a, 6], the covariances of the moments can always be expressed in terms of a set of centred 
covariances, i.e. corresponding to a limiting density with support centred at the origin. This allows us to write an explicit 
self-contained formula (28) for the limiting covariances of any one-cut /3-ensemble. Since the proof is based on the Dyson 
Coulomb gas analogy it is independent of the particular matrix realization of the ensemble and as such is simpler and more 
transparent than methods previously available. 

The classical /3-Gaussian, /3-Wishart and /3-Jacobi ensembles can be treated in a unified way. In this way we recover results 
previously established [24, 25] with proofs based on sparse matrix realization of /3-ensembles. In addition we have proved new 
formulae for the /3-Wishart and the /3-Jacobi ensembles. We mention that the issues of second order cumulants and joint behavior 
of spectral linear statistics have recently attracted attention in the context of the theory of quantum transport in mesoscopic 
cavities [19, 20, 30] as well as in practical spectral inference [18, 23]. 

We conclude by mentioning several related open problems. The challenge posed by higher cumulants is still open and the 
promising features of the 2D Coulomb gas approach suggest that the same line of reasoning could be pursued to investigate 
the universality of the o-point correlation kernel for v > 2. Concerning the first and second order cumulants, it is natural to 
look for the next to leading term corrections. In [35], Johansson showed how to compute the subleading terms for the averages 
at any /3 > 0 (for V polynomial). We are not aware of similar results for the o(N~ 2 ) corrections of the covariance structures. 
Coming back to leading order effects, it would be interesting to obtain an analogue of (25) for multi-cut /3-ensembles where the 
universality of the two-point smoothed kernel is more delicate [1, 2, 9], Asa last remark, we point out that the Coulomb gas 
approach is tailored to invariant ensembles. Wigner matrices (other than Gaussian) escape from the unitary invariance and this is 
a serious problem when discussing the n-point correlation kernels. For Wigner matrices it is known that, while the density of 
states is universal (the Wigner law), the higher-order correlation kernels are somewhat less universal [21] (they depend on more 
details of the entries distribution). The solution to these problems may well uncover yet another layer of “universality” features of 
“classical” RMT, a field from which many lessons are surely still to be learnt. 
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Appendix A: Gaussian Matrix integrals, Planar diagrams, and 2D Coulomb gas at /? = 2 

Enumeration of maps, i.e. graphs drawn on certain surfaces, random matrices and 2D Coulomb gases are inter-related topics in 
mathematical physics. The electrostatic analogy suggested by Wigner [56] and Dyson [27] in the late 1950’s has proved extremely 
useful in many branches of RMT. On the other hand, extensive elaborations on the connections between random matrices and 
the enumeration of maps did not begin until the works in the 1970’s by’t Hooft’s [48], Koplin et al. [37] and Brezin et al. [12] 
on the planar approximation in various field theories. Once this connection was brought to the surface, matrix models and field 
theoretical techniques have been exploited to make progress on unsolved graphical enumeration problems [8, 33, 34, 42]. (A 
comprehensive account can be found in the classical references [22, 57].) 

Diagrammatic techniques have been successfully applied to a few specific but important ensembles, namely the Gaussian 
ensembles, see for instance [53]. The crucial reason is that computing expectation values in this ensemble is tantamount to 
applying Wick’s formula in various settings; furthermore, when pairing the Gaussian entries, a nice geometrical structure emerges. 
These patterns have been rediscovered many times by different communities and mastering the relevant vast literature can be 
a daunting task for the uninitiated. The aim of this Appendix is twofold. Firstly, we summarize the general ideas connecting 
matrix integrals and planar enumeration of maps, and at the same time we try to clarify different technical terminologies from 
various communities. Secondly, we collect results available in the mathematical physics literature and compare them with our 
findings. Throughout this Appendix, we will concentrate on the GUE ensemble and on the diagrammatic techniques associated to 
the computation of its moments. For GUE, the joint law of the eigenvalues (1) is the canonical distribution of a 2D Coulomb gas 
confined in a harmonic potential V ( x ) = x 2 /A at inverse temperature /3 = 2. 
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ki = 4 , k2 -5 , k3 = 3 




FIG. 3: Left: An example of the mapping between NCPP and planar diagrams with a vertex (or planar maps with a star). In this example, the 
star has valence re = 8. This diagram arises in a i/? 8 -theory. In RMT this diagram contributes to the expectation (Gjv,s)- Right: This figure 
should clarify what happens when several moments are considered. We have a non-crossing connected pairing of three circles with a different 
number of points on each circle (note that the total number of points is even). The corresponding planar diagram has three vertices (v = 3) 
with different valences. It arises in a mixed field theory with cubic, quartic and quintic interactions. Diagrams of this type are counted when 
computing the third-order mixed cumulant Cs(Gn, 4 , Gat, 5 , Gn, 3 ) of the GUE ensemble. 

The idea that the large N asymptotic expansion of a matrix integral around a saddle point can be represented by diagrams 
identifiable with certain maps was pioneered by Brezin et al. [12]. In order to study the distribution of Gn,k, = TrS)^, for fixed re, 
it is customary to consider the function 

exp {-A^-Eat^s)} = (exp {-N 2 s G NjK }) , (Al) 

where hereafter the angle brackets denote the average with respect to the GUE measure. From the statistical independence of the 
Gaussian entries of Sat, it is known [12, 33] that Eat iK (s) admits a l/iV 2 -expansion 

E N>K (s) = Y / ^ g E i 9) (s), (A2) 

g> 0 

where, for all g > 0, Ei 9 ' 1 (s) is analytic in a neighbourhood of s = 0, 

V>1 

In the above expansion, Dif' 1 (v) is the number of inequivalent connected diagrams of genus g with v vertices of valence re. 
Therefore, in the large N limit only the g = 0 contribution survives, 

lim E n , k (s ) = £(°)( S ) = - E T^ 0) («)- (A4) 

N—too z ' v\ 

V>1 

The quantity E^\s) is often referred to as cumulant generating function of the re-th moment Ga/>. The coefficients in this 

expansion count the connected diagrams of a <p K -theory in which all non planar diagrams are omitted: ul!’ 1 (v) is the number of 
connected diagrams with v vertices all of valence re that can be drawn on a plane (genus <7 = 0) 6 . The probabilistic meaning of 
these numbers is clear from (Al): if we denote by C v (Gn,k, ■ ■ ■, Gn,k) the v-th cumulant of Gn,k., then by differentiating (Al) 
we arrive at 

lim N 2 ^C v (G NtK ,G N K ) = dM(v). (A5) 

N—>00 \ _ _/ 

v times 

For instance, the first two cumulants (average and variance) of Gn, k are given by the number of planar diagrams with one and 
two re-valent vertices, respectively, 

(Gat,«) = -0^(1) + o(l), (A6) 

iV 2 Var(G iv ,re) = D^\ 2) + o(l). (A7) 


6 


In the mathematical literature, these planar diagrams are also known as planar maps with v stars (each star has k, points). 
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We emphasize the origin of the two main features of these diagrams: the connectedness is intrinsic in the definition of E^ tK (s) 
(the logarithm of an exponential generating function); the planar approximation g = 0 emerges in the large N limit. Correction 
terms are encoded in E^ 91 (s) for higher genera g > 1, that is by considering k- vale lit connected diagrams l) 1 ,! 11 ( v) on higher 
genera surfaces (for instance the first correction comes from legal diagrams D^\v) on the torus, g = 1). In the one-vertex case 
(D^ (1) with k > 1), these planar diagrams are in bijection with non-crossing pair partitions (NCPP) of k points. A NCPP 
of k points can be represented by placing the numbers 1..... «: around a single circle and connecting numbers in pairs without 
crossings. Hence 


.0^(1) = =#= {non-crossing pairings of n points on a circles} . 

The theory of non-crossing partitions is an important topic in modern combinatorics and their role in random matrices has been 
elucidated by free probability [54, 55], with the notion of free cumulants [41, 45], 

One can extend the above identification to diagrams with more vertices v > 1 and write 

( v ) = ^{Non-crossing pairings of v circles with n points on each circle such that all circles are connected.} (A8) 


(See Fig. 3.) The connectedness requirement is clear and the non-crossing condition is the translation of the planarity condition 
for diagrams. Eqs. (A6) and (A7) thus connect average and variance of a single linear statistics Gjv, k of the Gaussian ensemble at 
(3 = 2 to the problem of non-crossing pairings on one circle, or on two circles, respectively. 

One can generalize the above arguments and consider the joint distribution of (Gm- jKi ,..., Gn,k„ ) f° r n >2 with distinct nfs 
in general. Using the notation k = («q, ..., K n ) and s = (s i,..., s n ), we may introduce the multidimensional analog of (Al) 

n 

exp {—N 2 En^(s)} = (exp {—N 2 SiG N , Ki }). (A9) 

i =1 

The large N limit of En^(s) will be the joint cumulant generating function of (Gn, Ki ,..., Gn, K ti ), 


Ef(s) = - 


E 


E 

v>l Vh...,v n >0 \ji =1 
uiH-1 -v n =v 


n 


(-*)* 

Vi'. 


D z\vV n ) 


(A10) 


where D^ r> (tq,..., v n ) is now the number of connected planar diagrams with iq Wi-valent’ vertices, V 2 ‘fi^-valent’ vertices, and 
so on. The relation with the theory of non-crossing partitions is now 

D%(vi,..., v n ) = ^{Non-crossing pairings of v = v± + ■ ■ ■ + v n circles with points on ui circles, 

vA.11 ) 

k -2 points on v-j circles, points on v n circles such that all circles are connected.} 


Clearly, no such pairing is possible if the total number of points (C” v t k, is odd. These numbers provide the leading order of 
the joint mixed cumulants (Gat jKi ,..., Gjv, k „ ) of order v = iq + • • • + v n , 


lim 

N—>oo 


N 2 ^C v (G n>Ki ,...,G 


N,ki y •• 


> G 


N,K n 


■,G N:K J = D { y\vv n ) . 


v\ times 


v n times 


(A12) 


From (A12), one can see that for the GUE ensemble the family {Gtv,k}k>o is asymptotically Gaussian (the case « = 0 is 
degenerate) as N goes to infinity. In the language of free probability, this means that the GUE ensemble has a second order limit 
distribution [16, 39, 40], i.e. the limits 

al= lim (Gjv.k) and {l/2)a s t = lim N 2 [(G N K G N e) - (G N K ) (G W f)] (A13) 

/V—Vrv) 5 /V—Vrvi 


exist for all k, ( > 0, and the higher-order cumulants decay faster to zero for large N. 

Brezin et al. [12] considered the planar approximation for n = 3 and k = 4 (quartic and cubic vertices). They computed exactly 
the cumulant generating functions E^ (s) and E^ ( s ) thus providing the. full family of cumulant of TrS 3 and TrS 4 to leading 
order in N. For the quartic interaction, the higher genus corrections E^\s) and E^\s) have been computed explicitly in [8]. 
The problem of mixed interaction (diagrams whose vertices have mixed valence) has been less explored even in the planar regime 
g = 0 — the case of the combined cubic and quartic interactions, i.e. the computation of E^ J j (,sq, s 2 ) was proposed in [12] but 
was not worked out explicitly. Later, Harer and Zagier [31] computed the exact finite N value of (Gw K ) (i-e- the linear term in s 
of Ejsr tKj (s) in (Al) for all n). 
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FIG. 4: Non-crossing pairings of two circles with k points on the first circle and l points on the second circle, such that the two circles are 
connected. The number of such pairings is equal to D} 0 )(1, 1) = (l/2)a® t , i.e. limjv-Kxj Cov(TrSjv, TrS^r) with Sjv belonging to the GUE. 
We have (1/2)014 = T ( 1 / 2 ) 02,2 = 2, (1/2)04,3 = 3, (1/2)034 = 12 according to Eq. (42) and Table I. 


Our results (when restricted to the GUE ensemble) provide a complete answer to the following problem: for any finite family 
of power traces (Gjv >Ki ,..., Gjv, Kn ) compute the joint cumulant generating function E l 0 '* (s) in ( A 10) up to quadratic terms, 


4 0) (^) = -E E 


l<v<2 vi ,...,v n >0 

■uiH- \~Vn=V 


n 

.i=l 


(-Si ) 1 


■■■,V n ) + ■ 


(A14) 


= E aE a. i) + 


(A15) 


i,j =1 


This corresponds to considering the contribution from all types of planar diagrams (every kind of combined interactions) with 
at most two vertices (of arbitrary valence). Here (1) and D^ Kj (1,1) give the leading order behaviour of the averages and 
covariances of the family (Gjv) Kl , ..., Gn,k„)- The results of Section IV A, restricted to /3 = 2, give 

lim ( G N , K )=Dl°\l) = al 

N —100 

lim N 2 Cov(G n K , G N i) = D[ 0) e ( 1,1) = (1/2 )o% t 

The explicit formulae were given in (40) and (42), respectively. Therefore, from the above discussion we have 

a K, — #{Non-crossing pairings of k points on the circle}, 

(1/2)a® i = #{Non-crossing pairings of two circles with n points on the 1st circle 
and i points on the 2nd circle such that the two circles are connected} 

In Fig. 4 we report few such diagrams. 

The difference of our method compared to previous works on the planar approximation in field theories is the following. 
Instead of expanding the cumulant generating function of a fixed Gjv, k (or a fixed rc-uple (G \;. K] ,..., C!\;. Kn )) in the number 
of vertices v as in (A4)-(A10). in this paper we have considered the whole family of generating functions of a generic n-uple 
(Gjv Kl , • ■ ■, Gjv ;Kri ) up to quadratic terms. In diagrammatic terms, instead of an expansion of planar diagrams with fixed valence 
k in the number of vertices v, we would obtain a series in the valences of vertices n and £ of planar diagrams with at most two 
vertices v < 2. 

We mention that the combinatorial problem associated to a% and (1/2)a® t was solved by Tutte [52] as a particular case of a 
more general problem on the number of so-called slicings of a band. We note that Tutte’s formula enumerates general slicings, 
with a strong restriction on the parity of the vertices. In the language of RMT, Tutte provided a formula for the generic mixed 
cumulant D^\v 1 ,... ,v n ) of (Gat jKi ,..., Gn,k „) (to leading order in N ), with the condition that Ki,..., kn should all be 
even 7 . We find interesting to note that this obstruction (that might appear unjustified in the combinatorial problem) can be 
explained in field-theory with the instability of any tp K -interaction for k odd. 


(A18) 

(A19) 


(Alb) 

(A17) 


7 In the case of only two bounding curves (our case), Tutte managed to obtain a formula that comprises the cases tt and £ both even or both odd. No orientation is 
involved in his problem, hence in his formula there is a factor of 1/2. 
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Appendix B: Outline of the proof of (42) 

It is convenient to cast or K f in (42) for re,/: > 0 in term of Euler Gamma functions 

16T(re)r(f) 


vS = J (K + *)r(f) r(f) 
reM^-xrd) 2 r(|) 

7T 2 4 (re + /')r(re)r(f) 


for £ and re even 


for i and re, odd . 


(Bl) 


Setting re = 2 r and 1 = 2s (in the even-even case) and re = 2r + 1 and £ = 2s + 1 (in the odd-odd case), we can evaluate the 
double sum in (25) as 


E = E 


8T(2r)T(2s) 


r c 2s + E 


(2r + l)(2s + l)4 2r+2s+1 r (yti) 2 r (^±±) 2 
7r 2 (2r + 2s + 2)r(2r + l)T(2s + 1) 


2r+l a2s+1 


k,/= 1 r,s=l ( r + S 6 7 8 9 10 ) r ( r ) r ( S ) r,s=0 

The two series on the r.h.s. can be computed in closed form using the following identities 

i r°° 

d te~ t(r+s) , 


r + s 


V F ( 2r ) r„ c -t/2^2r = e tz 2 

“t r (r) 2 “ (e-* (e* - 4 2 2 )) 3/2 ’ 

V ( 2r + 1 ) r ( g E 1 ) (/|~r-«) 2r =_ — _ 

^ r(2r + 1) (e 2t - 4 z 2 ) y/\ - 4e" 2 ^ 2 


Therefore we can rewrite the double sum in integral form 


E <t zK c l = s / df 


e-‘C 2 


AC,^=1 


o (e _t (e* — 4^ 2 )) 3 ^ 2 (e _t (e* — 4^ 2 )) 3 ^ 2 


+ ^C/ 


7re 


7re 


( e 2 * _ 4^2) yj \ _ 4 e 2 * z 2 ( e 2t _ 4^2 ^ — 4e _2t </ 2 

The integrals can be carried out in closed form, and after simplifications we obtain 

C,z (y 1 - 4CV1 ~ 40 2 + 4Cx - 1 ) 


E = - 


K ,, =1 - o 2 

which (after multiplying by 1//3) coincides with F± 2 {z, Q as given in (41). 


for -1/2 <z,C< 1/2 


c 2s+i . (B2) 


(B3) 

(B4) 

(B5) 


(B6) 


(B7) 
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